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Abstract 

The methods of conformal field theory are used to obtain the series of exact solu- 
tions of the fundamental equations of the theory of turbulence. The basic conjecture, 
proved to be self-consistent ,is the conformal invariance of the inertial range. The 
resulting physical picture is different from the standard one , since the enstrophy 
transfer is catalyzed by the large scale motions. The theory gives some unambiguous 
predictions for the correlations in the inertial range. 



1 Introduction 



During the last thirty years we witnessed an amazing unification of physical ideas. Con- 
cepts and notions of the totally different regions of physics appeared to be almost isomor- 
phic. Spontaneous symmetry breaking and renormalization group are the most notable 
examples. 

In this paper I extend slightly this set of isomorphic ideas. Namely, I apply methods 
of conformal field theory ( CFT ) to the problem of two-dimensional turbulence. Short 
version of this paper has been published earlier [1]. 

The major puzzle in the theory of turbulence is the following. It is commonly believed, 
that at large Reynolds numbers we are dealing with the stationary statistical regime, 
describing velocity distribution. In other words, there should exist a time-independent 
probability P = P[ v a (x) ] (where v a (x) ) is the velocity). This probability should 
commute with the equations of motion, or: 

/*G^))*" (l) = (1) 

(where it is understood that we express v a through v a by the use of Navier - Stokes 
equations). At large Reynolds numbers viscosity can be neglected and hydrodynamics 
becomes a hamiltonian theory with the hamiltonian: 

H = J dx\vl (2) 

If so, the equation (1) means that P is the integral of motion for this system. Since it 
is believed that in general there are no non-trivial integrals for the system (2) (in two- 
dimensions there are, but these are irrelevant for our discussion) we could erroneously 
conclude that: 

P =P(H) = exp{-(3H} (3) 

(the last equality follows from the usual additivity assumption). This is the Gibbs dis- 
tribution for the temperature and it clearly does not describe the turbulent flow, in 
which we expect permanent energy flux. 

The resolution of this puzzle is that in fact there are extra integrals of motion which 
form P[ v a ]. They are, however, highly non-local and non-polynomial. It would be 
difficult and unnecessary to write them down explicitly, because there still be a problem 
of averaging with respect to very complicated distributions. In what follows we will 
mostly use the analogue of the equation (1) applied directly to the correlation functions. 
Nevertheless we will briefly comment on the origin of the extra integral of motions for the 
systems with weak coupling. In this case they are just naive field-theoretic constants of 
motion formed of " in" and "out" operators. When expressed in terms of the Heisenberg 
fields they become non-local, and describe the distributions of weak turbulence. 
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The central point of this paper is, however, quite different. We will consider the 
equations for the correlation functions: 



< v ai (x ± ) v a2 (x 2 ) • • • > + < v ai (xi)v a2 (x 2 ) - • • >H = 



(4) 



(where v a is expressed in terms of v by the use of equations of motion). These are standard 
equations (The Hopf equations [2] ) which express N - point functions in terms of the N 
+ 1 - point function. The existing procedures for dealing with these equations are based 
on some kind of closure hypothesis expressing, say the 4 - point function as a square of 2 
- point function. 

We take a different approach. Namely, we try to satisfy (4) exactly, by the use of 
conformal field theory, assuming that the developed turbulence in the inertial range pos- 
sesses infinite conformal symmetry. There is no proof that this must necessarily the case, 
and that our conformal solutions are the only possible ones. What we have done is only 
a successful ansatz which solves exactly eqs (4). 



2 Chains of equations solved by conformal field the- 
ories. 



In this section I will sum up some basic properties of CFT, needed below as a represen- 
tative example, I will use the most familiar case of the Ising model. The Ising model is 
a tp A -field theory [3]. That means that if we treat (p as a fluctuating variable, it satisfies 
the equation: 



implying the relations between N - point functions of <p and N + 3 - point functions, 
known as Schwinger - Dyson equations. The constant fio must be fine - tuned in order to 
put the system to the critical point, or in order to allow (5) to have conformal invariant 
solutions. Now let us try to answer the basic question of this section : how CFT solves eqs 
(5)? In CFT we are dealing directly with the sets of correlation functions at the critical 
point. The basic assumptions [4] are that there exists a set of the so called primary 
operators {Ok(x)} with anomalous dimensions (^Ak,AkJ (referring to the rescaling of 
variables z = X\ + ix 2 and z = x 1 — ix 2 . It is assumed that under arbitrary analytic 
transformation z =>- f (z) ,Ok transform as: 



This transformation is generated by the analytic energy- momentum tensors T(z) and T(z) 
(related to the tensor T a p (x) as T, T = Tn — T 22 ± 2iT\ 2 ). Conformal invariance implies 
that the trace of the energy-momentum tensor is zero. 




(5) 
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The basic result of [4] was that together with primary operators, CFT is bound to 
contain the so-called secondary or descendent operators. They come from the conformal 
deformations of the primary ones. In order to describe them it is convenient to introduce 
the Virasoro algebra: 

L n = j dzz n+1 T(z) (7) 
[L n ,L m ] = (n - m) L n+m + ^-n (n 2 - l) 5 n+mfi (8) 



k 



The secondary operators are given by: O^'" u 1 "' n = L_ ni . . . L_ n; L_„ tl . . . L_ mn O, 
The correlation functions of the secondaries are expressed by the differential relations 
through the correlators of the primaries [4]. 

The second basic feature of CFT is the operator product expansion. It says that 

O k (z)O l (0) = ]T / fe 7(^) Am ~ Afc ~ A; O m (0) + secondaries (9) 

The structure constants satisfy simple (in principal) consistency relation coming from the 
associativity of the algebra. Thus, the classification problem of CFT is reminiscent of the 
classification problem of the Lie algebras, but is vastly more complicated. 

However, in [4] we have found a simplest set of CFT's, the so called minimal models. 
In these models the number of the primaries is finite (but the total number of operators 
is, of course, infinite). The Ising model is the minimal model of the type (3,4) with only 
3 primary operators. They are unit operator I , spin a and energy density e. The fusion 
rules (9) in this case are: 

[a] [a] = [I] + [e] (10) 
[a] [e] = [a] (11) 

m = [j] (12) 

The brackets in these symbolic formulas mean the " conformal class" , namely the primary 
operator itself together with all it descendents. 

Now we are coming to the central point. Let us identify the field ip of the p^ theory 
with the field a of the conformal field theory. We have to examine the quantities 
< p 3 (x) <p (yi) • • • tp (y n ) >. However, the conformal field theory is applicable only when 
all points are well separated compared to the lattice spacing. The way out of this problem 
is to introduce the point - splitted definition of the ip 3 (x). Let us examine the object: 

TETa ( x + l) a ( x -l) a ( x ) = " a3 ( x ) " = $ ( x ' a ) ( 13 ) 



The symbol lim here means that we first of all average over directions of a and then take 
\a\ to be much smaller then all other distances \x — yi\ and \yi — yj\. Keeping it at first 
much larger than the lattice spacing, we can apply the fusion rules (9). Since 

[a] [a] [a] = [a] (14) 



3 



We obtain: 

$(x,a) « \a\~ iA °\ Cl a(x) + c 2 a 2 d 2 o(x) + . . .] (15) 

Although this equation is correct only for \a\ » I (where 1 is a lattice spacing), with 
any natural definition of the theory it can be extrapolated to \a\ ~ /. The reason is that 
with any isotropic regularization at the distance we will get instead of (15): 

$(x,a)=r 4A °f(^pja(x) + ... (16) 

(where precise form of f depends on the way we regulate the theory). As a result, if we take 
a < I we still have the equation (15), and hence the basic field equation (5) is satisfied. 
All this means one simple thing -field equations of motion, when read from the right to the 
left are nothing but operator product expansions. This fact was well understood from the 
very beginning of the field theory for critical phenomena. Point-splitting and directional 
averaging over a is also nothing new - already in quantum electro-dynamics it has been 
necessary to introduce the isotropic cut-off in the momentum space, which is equivalent 
to the procedure we discuss. 

To summarize - in this section we have seen how the structure of the operator product 
expansion reflects equations, satisfied by the correlation functions. 



3 Solutions of the inviscid Hopf equations by the con- 
formal field theory. 

Let us turn now to the Hopf equations (4) and try to solve it in the spirit of the proceeding 
section. In two dimensions it is convenient to introduce vorticity u and the stream function 
■0 given by 

v a (x) = e a(3 dpij) lo(x) = e a pd a vp = d 2 ip (17) 

They satisfy Navier - Stokes equations: 

Co + e a pd a ip dpd 2 il) = vd 2 uo (18) 

If the stirring force is present, it must be added to the RHS of (18). If we assume, that for 
the large Reynolds numbers there exist the inertial range of scales in which both viscosity 
and the stirring force are negligible, we can as a first step examine the inviscid Hopf 
equation: 

<iv(x 1 ) u(x 2 )---> + <u(x 1 ) lj(x 2 ) ••■> + ••• =0 (19) 

w(x) = - "e^d^dpd 2 ^ (20) 

Here, just as in the previous section, we must be careful in defining the correlators at 
the coinciding points. The reason there have been tied to the existence of the lattice 
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cut-off (where conformal theory fails). The reason here is analogous - we expect that in 
the momentum space only a; (A;) with < \k\ < - should satisfy (19). Here a is the 
ultraviolet cut-off proportional, as we see later, to some power of the viscosity v while L 
is the infrared cut-off, dictated by the scale of the stirring force. The UV - cut - off in the 
momentum space means that we must use the point- splitted definition of the product in 
(20) __ 

"e a(3 d a i; (x) d p d 2 i) (x) " = lim e a t,d a 1> (x + a) dpd 2 ^ (x) (21) 

Let us assume now, that ip is a primary operator of some yet unknown conformal field 
theory and compute the RHS of (21) by the use of operator algebra. Suppose that we 
have the structure: 

ip (x + a) ip (x) = (aa) ^ {<p(x) + descendents} (22) 

Here we introduced the operator <fi which is the minimal dimension operator in the op- 
erator product (22). In unitary theories (as in the Ising model) it is a unit operator. 
However, only Gibbs states (from which " nothing disappears" ) are described by the uni- 
tary theories. In turbulence we have "flux states" from which conserved quantities leak 
away, and it is natural to expect that they will be described by the non-unitary theories. 
In this case most of operators have negative dimensions and (j) is non-trivial. It is clear 
that when the operation (21) is applied to (22), the leading term in the expansion gives 
zero. This happens because the result must be pseudoscalar and we can't form it out of 
(p and its derivatives. Hence we have to consider subleading descendents in (22). They 
have the general structure: 

descendents = c {n}{m} L_ ni . . . L^ nk L- mi . . . L_ mi (f) (x) 

(aa) A *~ 2A ^a£ n a£ m (23) 

Differentiation and directional averaging in (21) the leading term in (23): 

"e a pd a ip (x) d p d 2 i) (x) " = const (aa) A ^ 2A ^ [L_ 2 l\ - L^L 2 ^} 4> (24) 

This one of the basic formulas of this work. The reason, why this descendent appeared 
in (24) is quite simple. The LHS of (24) is pseudoscalar and hence it changes sign under 
complex conjugation. The lowest descendent which has the same property is the one we 
find in (24). In deriving (24) we accounted for the UV - divergency, but disregarded the 
IR one. We will return to the infrared side of the problem below. 

Let us assume for the moment that IR - divergencies are not present. If so, we derive 
a rather amazing conclusion from the eq.(24). Namely, there are two possibilities. The 
first one is A^ < 2 A,/, (negative:defect of dimensions"). In this case the only way to 
satisfy (19) would be to assume that the operator 

n = (L 2 L\ - L^L 2 ^) cj) (25) 
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is either zero or is a symmetry of the underlying CFT. The latter statement just stresses 
the fact that the relation (19) implies that the infinitesimal variation of uo , Co = Q 
must not change the correlation functions, or as in the case of decaying turbulence rescale 
them. It is quite possible that non-trivial examples of such symmetries exist. However 
if we restrict ourselves with the simple minimal models, we will have to conclude that 
the RHS of (25) must be zero. This is easy to achieve if we require that the minimal 
dimension operator is degenerate on the level two. That means that it satisfies the 
equation [4]: 

( L - - wi^rf-) * = (26) 

If so, n — . The simplest example, in which this relation is satisfied is the minimal 
model (2 , 5) with = A^ = — | 

Even more surprising is the fact, that if A^ > 2A^ then no extra restrictions are 
needed because as a — > the RHS of (25) is zero, and the inviscid Hoft equation is 
satisfied. So, any CFT with the positive defect of dimensions solves the Hopf equations. 

This is a "turbulent" counter part of the fact that there are continuosly many parasite 
solutions of the static Euler equations. Let us remember that if we look at the " laminary" 
solutions of 

e af} d a ^d d 2 il) = (27) 

we find that if d 2 i(j = f (tp) with arbitrary / (ip) then (27) is satisfied. It is well known, 
however, that matching with the viscous region by the boundary layer consideration re- 
moves this ambiguity. This analogy makes it clear that the crucial point for selecting 
correct solutions of the inviscid Hopf equations must be matching relations at the UV 
region of the cut-off momenta, where viscosity takes over. It is also clear, that since we 
are dealing with negative dimensions, we have to analyze the infrared diver gencies.We 
can say that in turbulence there are two boundary layers in the " momentum space" , each 
requiring a separate theory. In the following sections we will discuss some partial results 
in this direction. 



4 The infrared problem 

The CFT's which are most interesting for the turbulence problem contain negative di- 
mensions. That means that the formal correlation functions are positive powers of the 
distance. It is clear that we have to be very careful in the infrared region and treat it 
separately. Let us begin with the 2 - point function < (0) (x) >. When the dimension 
A^ < 0, CFT gives 

(conf ) . , . | 

<il>(0)if>(x) > = -N 4|A *! (28) 

This is clearly unphysical. The correct prescription should take into account the fact, that 
only momenta in the inertial range are described by CFT. So, we should expect that, if 
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we take the Fourier transform of (28) 

<ip{k) V> (-k) > = const | fc | 2 + 4 | A ^| (29) 

this result must be used only if L _1 < |fc| < a -1 . We do not know the contribution from 
the |fc| ~ L -1 . The only thing which can be said is that when returning to the x space it 
should give expressions, analytic in x 2 . So, we have: 

< V (0) i> (x) >~ J > l i e ikx \k\- 2 - 4 ^ ~ ( Cl L 4|A ^ + ...)- |x| 4|A ^ (30) 

with the first bracket representing the contribution of the infrared modes. When we turn 
to the multi - point functions, the situation is essentially the same. When we examine 
their momentum space form: 

G (fci • • • k N ) = <i)(k 1 )---i)(k N )> (Ekj — 0) (31) 

we will assume that the CFT formulas for this quantity work,provided that: 

L- 1 < \ki\ < a' 1 (32) 

and 

L' 1 < \k a + ■ ■ ■ + k u \ < a' 1 (33) 

The (33) condition means that we have no small momenta transferred in any channel. 
Again, in the coordinate space all that means that the physical N - point function is 
equal to the conformal one plus the terms, analytic in some (x, — Xj). Notice, that if we 
understand the momentum integrals as analytically continued in (A^,) from the region, 
where they converge, we get precisely conformal answer in the x - space. In order to get the 
physical correlator within this analytic regularization, one has to add 5 - functions in the 
k - space, representing the infrared modes. For instance, within the analytic prescription, 
we have instead of (29) 

<V(*0V(-*)> = Jj^^i + a lL ^5(k)+a 2 L^- 2 ^ 2 S(k) + ••• (34) 

When Fourier - transformed, (34) gives again (30). 

So, the convenient way to summarize this situation is to say, that we have some sort 
of Bose - condensate in the momentum space, formed by the large scale motions. The 
physical correlators differ from the conformal ones by the condensate terms (we will call 
them S- terms later). The 5 - terms occur whenever any sum of momenta in the correlator 
is equal to zero. In the coordinate space they are represented by the terms analytic in 
(xi — Xj). Occurrence of the 5- terms is not surprising , because we have strirring forces 
acting at k=0. The main problem with them is that we do not know yet how to determine 
their form from the dynamics. The other related problem, which we will consider now, is 
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whether they destroy our solution of the Hopf equations. The danger comes from the fact 
that in the sect. 3 we dealt with the purely conformal correlators, and this deri vation has 
to be reconsidered. In order to do that, let us rewrite our basic equation in the momentum 
space: 

<co(q)co(f 1 )---co(f n )> = -jd 2 k[k,q]((q-k) 2 -k 2 ) < if) (k) if) (q - k) co (A) • • • u (/„) > 

(35) 

Generally speaking, this integral has infrared divergency as k — > (or k — > q ). The 
results of the proceeding section guarantee that if we treat this IR divergency by the 
analytic regularization, the RHS of (35) will be zero, as before. 

The problem now is to study what happens when we simply impose the IR cut off in 
(35). In order to do this, we need to know the behavior of the correlations functions when 
one of the momenta is much smaller then the others. Let us first evaluate the divergencies 
at k — > 0. This limit in the coordinate space is dominated by the large R limit in the 
Fourier transform: 

< if)(k) if){q-k)co{h)---co{f N ) > = Jd 2 Re lkR <if)(R)if)(0)u)(f 1 )---u(f N )> 

(36) 

The large R behavior is governed by the fusion rule, where we replace all operators except 
if) (R) by if) (0) . Using notations of [4] , we have the formula: 

<if)(R)i)(0)co(f 1 )...> = ]T< 0|V>OR) |A + iV >< A + N\if)(0)co(f 1 )...\0> 

N 

(37) 

where the states |A + iV > correspond to the secondary operators of the if). Using the 
fact that 

< 0\ 1 f J (R)\A + N>oo l ^ w (38) 
we obtain the following leading terms in the asymptotic expansion 

< ip (k) if) (q — k) co (/i) ■■■co (f N ) < (q) cofi • • • w (f N ) |0 > 

+ c 2 ^j^<A\L 1 if)(q)u(f l )---u(f N ) 0>+c. C J + ..- (39) 

When we substitute (39) into (35), the first leading term drops out after integration over 
directions of k the second term does contribute. Simple computation gives: 

<d)(q) co (/i) ■■-co (f n ) >= constL 2+x {q + < A\L ± co (q)---co (/„) | > -c.c.} + O (L A ) 

(40) 

This leading term can be interpreted as fluctuating transport by the large scale eddies. 
Indeed, the relation (40) is equivalent to the operator equation: 

co (x) = L 2+x u a d a co (x) (41) 
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Here we introduced a new x - independent operators u±hy the rule: 

< 0\u + u (xi) . . . u (x n ) |0 > = % < A^\LiUJ (xi) . . .uj (x n ) |0 > (42) 

andw_ = u* + . It is clear, that although this infrared counterterm makes Co 7^ 0, the Hopf 
equation is still satisfied, due to the translation invariance. As we go to the order of 
L A counterterms we find more complicated forms. By the same technic it is possible to 
show, that the complete divergencies give the following equation: 

Co = u a d a ui + h af3 d a dpi) (43) 

with the 

< 0\h ++ u (xi) •• -u (x N ) |0 > ~ < A| (l 2 + XL^\ u (xi) ■ ■ -u (x N ) |0 > (44) 

(where h a p- is a traceless tensor). This contribution does not cancel in the Hopf equation. 
There are also other infrared counterterm with this property, occurring from the region 
in the (45) where: \k + fn + • • • + fu\ — > (for some choice of {fi} ). All that means 
that if we had defined the physical correlation functions by simply cutting of the momen- 
tum integrals in the Fourier transform, the Hopf equation would be destroyed. On the 
other hand, if we had defined all integrals by analytic regularization it would mean the 
identification of the physical and conformal correlators. In this case the Hopf equations 
is satisfied, but the positivity is broken. The way out of this problem is to find the 5 - 
terms which do not spoil the Hopf equation but restore the positivity. The most blatant 
violation of positivity is easily cured, if we postulate that any correlation function contains 
the constant part, coming from the zero modes: 

< O nl (x 1 ) O n2 (x 2 ) • • • CW (x N ) > phys = c nl ■ ■ ■ n N (L~ 2 ^ A -) + < O nl ■ ■ ■ O nN > conf 

(45) 

(where C - are some constants). This type of 5- term does not contribute to the Hopf 
equations at non - zero momenta, and at the same time prevails over the negative con- 
formal contributions, since \xi — Xj\ L. That doesn't mean that there could be no 
other IR counterterms. In fact it is possible to show that there are many types of the 
5 - terms, consistent with the Hopf equations at g.; 7^ 0. In order to determine them we 
need to include the stirring force terms (which arise at qi = 0) and to use the matching 
relations. This is the problem for the future work. At present it is not quite clear to what 
extent the 5- terms are universal i.e. independent of the large scale structures. 

There is another related issue which we should discuss in this section.^ It is the 
question of the vacuum expectation values of different operators. When defining our 
conformal contributions, we have implicitly assumed that all these expectation values 
(VEV) are zero. In general, however we have: 

< O n (x) > = C n ■ L- 2A " (46) 

1 I am grateful to A. Zamolodchicov for drawing my attention to it 
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If A n < , and c n 7^ these VEV considerably modify correlation functions, as was 
noticed by Al. Zamolodchikov[5 ]. The reason for that is simple. From the OPE we get: 



< V (r) i) (0) > = r~ 4A ^ < / > + r 2 ( A ^ 2 M < >+••• = r 4 l A *l \ 1+ < > • 



The second term in (48) is the dominating one. The question, whether VEV are non-zero 
is determined by the infrared boundary conditions. Our assumption that c n = is not 
fundamental. Actually our consideration of the Hopf equations remains intact even if 
c n 0. The only thing which becomes more complicated is the infrared divergency. Also, 
the UV - matching of the next section requires modified consideration. 

The question of the c n values belongs to the same category of the IR problems which 
we left for the future work. At present the situation resembles quantum chromodynamics, 
where we have not solved the IR problem, but can, nevertheless, test its small distance 
behaviour. 

Let us summarize. The true physical correlators can be obtained from the conformal 
ones either by defining the momentum space integrals in the sense of analytic regular- 
ization or in the usual sense, with the IR cut - off. In both cases it is necessary to add 
to these expressions the 5-terms, containing 5 functions of certain external or transferred 
momenta. These terms in the coordinate space contain the analytic dependence on cer- 
tain distances. Precise form of the 8 - terms requires for its definition the matching with 
the stirring force - the task not accomplished in this paper. 



In the above considerations we accounted for viscosity in a somewhat symbolic sense. 
Namely, we have used the inviscid Hopf equations, but assumed that at the very high 
momenta our power-like correlation functions start to drop rapidly because of the viscosity, 
which was assumed as the origin of the ultra- violet cut-off. 

It is clear that the precise nature of this cut-off must be important. For instance, 
we can imagine an "elastic" cut-off which preserves the hamiltonian structure and the 
"inelastic" one which introduces dissipation. It is clear that they should correspond to 
different physics. 

In this section we' 11 try to impose the dissipation condition on the theory. The 
basic idea is roughly the same as in the classical Kolmogorov treatment of turbulence. 
Namely, we will assume (self-consistently) that integrals of motion are dissipated at the 
UV cut-off, and in the inertial range of momenta they are just transferred. It implies that 
the flux in the momentum space of conserved integrals must be constant. One can say 
that while Gibbs distributions are uniform on the surfaces of fixed values of conserved 
quantities, the turbulent distributions are located on the surfaces of the constant fluxes 




(47) 



5 The Flux States 
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of the corresponding quantities. 

The constant flux in the momentum space is also a familiar object in the other part 
of physics. It is responsible for the anomalies in the quantum field theory. When one 
considers, for example, massless fermions in the electromagnetic field the chirality (the 
number of left minus the number of right Dirac particles), is not conserved due to axial 
anomaly. This happens because the ultra-violet regularization breaks conservation of 
the axial current. After the chirality is injected at the UV cut-off, it propagates in the 
momentum space to the physical region. There is, therefore, a clear and useful analogy 
with what we are discussing. 

It has been noticed by Kraichnan [6], that in two dimensions the most important flux 
is that of enstrophy (we will comment on that later). Let us derive the constraint on 
the theory which follows from the enstrophy conservation. Consider first the enstrophy 
contained in the modes with definite momentum: 



<u(k) uj (-k) > (48 



h{k) = 

Without viscosity and external force, h(Jcj would satisfy the continuity equation in the 
momentum space. When these factors are accounted for, the equation becomes: 

h(k) + ^-4 h) = vk 2 h(k) + <S>(k) (49) 

OKi 

Here $ (Jtj is a contribution of the external forces which is non-zero only for \k\ ~ 

Let us notice also, that in isotropic turbulence the only non-zero component of is 
the radial one (k). Let us now chose the momentum q lying in the inertial range and 
integrate over \k\ > \q\. Since in the steady state h(k) =0 we get 

- jCO ( q ) = v [ k 2 h (k) d 2 k (50) 
J\k\>\ q \ 

Now comes an important point - suppose that the RHS of (51) is UV- divergent and 
defined by \k\ ~ ^. Then the q -dependence of RHS can be neglected. As a result, we 
obtain: 

J {h) (q) « -v J k 2 h (k)d 2 k = const (51) 

for jr <C \q\ <C i. The meaning of this constant flux condition is transparent. It just says 
that enstrophy is dissipated only at \q\ ~ ^ and hence, the conservation law implies that 
its flux must be constant in the inertial range. This is one of matching conditions we 
must impose on the inviscid solutions. If we express u (q) by means of eq. (35) we get: 

jW(q)= -[ < uj (k)u (-k) > = [ <u(k)u(-k)> (52) 

J\k\>q J\k\<q 

(we used here the fact, that (35) conserves total enstrophy). We see from here that all 
k — s from the inertial range give zero contribution. At the same time, using our general 
conjecture concerning the infrared terms (46) we obtain: 
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The right hand side of (51) is independent of L and thus we obtain the matching condition: 



A w + A* = (54) 

If we had used the first term in ( 52), where k ~ ^ are absent, the IR contribution 
would come from the integral (35) for Co , yielding the same result. If we recall, that: 
A w = A^ + 1 and A^ = A^ + 2 we get a following condition for the enstrophy flux 
state: 

M M = M + • • • 

A4, + A + 3 = (55) 

This result heavily depends on our conjectures about IR - contributions, namely on the 
fact that < ujuj > infrared part is non-zero. Had we considered the energy flux state, the 
formula ( 55) would be replaced by: 

A^ + A^ + 2 = (56) 

Apart from the energy and the enstrophy there are also higher conserved integrals of the 
type: 

I n = Ju n (x) d 2 x (57) 

Before commenting on these integrals, let us explain why the formulas (56) and (57) 
do not contradict each other. (After all, energy and enstrophy are both conserved in the 
inviscid system) . The apparent contradiction is removed by the fact that in the enstrophy 
flux state the integral for the energy dissipation, analogous to (50) is given by: 

J&(q)=v[ d 2 k <u{k)u (—k) > (58) 

J\k\>q 

Contrary to (50), (58) is not UV - divergent and the dissipation is scale - dependent in 
the inert ial range. 

In this case we can't use (51) anymore. No obvious matching condition arise in this 
case, since (58) tells us that as v — > the energy flux tends to zero, while the enstrophy 
flux persists. This is essentially the picture, advocated by Kraichnan [6]. 

As we turn to the integrals I n we have to decide, whether they have nonzero flux and 
whether this gives us any new information. We are unable at this point to give complete 
analyses of this question and will present some non-rigorous estimates, leaving complete 
resolution for the future work. We have: 



<V 2 cu(x)u n - 1 (x) >d 2 x (59) 

The main problem is how to regularize this expression. We must use OPE in order to 
define u) n ~ l .If n is even, which is imposed by parity conservation, we get 

uj n ~ l (x)ooV(x) + ■■■ (60) 
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and result: 

4 oc J < Au V > oc J (e) (61) 

We picked up if) operator only in(61) since due to orthogonality of primary operators, all 
others will not contribute. This implies that the I n - transfer is zero in the inviscid limit 
and there are no extra constraints on the theory coming from it (apart, may be from some 
inequalities). However, this line of arguments must be considered in more details (with 
careful point - splitted definition of all parties involved, before the last conclusion could 
be trusted). In particular,the other groupings of to -s may involve some other operators. 
Much work here is to be done. 

An interesting related question is that of the inverse cascades of higher integrals. Q 
Strong infrared divergency in their dissipation makes them very similar to the energy 
and opens this possibility for I n . Whether this is the case is an unsolved question. 
Unfortunately, any kind of the naive pseudophysical arguments is misleading here. 

To conclude this section, let us discuss the physical picture of the flux state. It is not 
exactly the Kolmogorov's one. We have constant flux of conserved quantity through the 
inertial range. In the Kolmogorov's mechanism the transfer occurs by the interaction of 
three modes with all wave vectors lying in the inertial range. Such transfer is absent in 
our picture. Instead we have the triad interaction in which one of the modes is infrared. 
We can say that the infrared (large scale) modes serve as catalyzers for the flux. This is 
an important qualitative prediction of our theory, which can be tested. 



6 Possible spectra of conformal turbulence. 

Let us discuss solutions of conformal turbulence with constant enstrophy flux. The analy- 
ses of the previous sections indicates that we must start from the CFT, satisfying certain 
constrains which follow from the Hopf equation and constant flux condition. The first 
constraint tells us that if we have the fusion rule: 

m w\ = M + ■■• (62) 

where is a minimal dimension operator in the product (62), then it must be either 
A^ > 2A^p or - operator must be degenerate at the level two. 

The second constraint depends on the assumption about vacuum expectation values. 
The simplest version of the theory, when these values are zero gives: 

A + A^ = -3 (63) 

More options arise if we assume that if < O n > ^ In this case the correlation func- 
tions toco and 6joj are determined not by the unit operator but by the minimal dimension 

2 1 acknowledge useful discussion of this point with G. Falcovich and A. Hanany 
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operators. As a result, if we look at OPE, we find: 

< CO (z) CO (0) >OC {zz) A *~ 2Aul <<p> 

A w = A^ + 1 (64) 

and 

< co co >oc (zz)- iA " +A ^ (zz) Ax (65) 

where x is a minimal dimension operator in the product [ip] [ip] [ip]. So, in this case the 
constant flux condition is replaced by: 

+ A^, - A x = -3 (66) 

while the energy spectrum is given by: 

E (k) ~ k AA *~ 2A * +1 (67) 

One can also consider intermediate possibilities when some of the operators have vacuum 
averages and some do not. 

At present we do not have clear idea how the choice should be made. It may depend 
on the external conditions and the way turbulence is generated. 

Below we will examine a simplest possible model, the "minimal" minimal model. By 
that we mean the theory which satisfies our requirements and has the smallest number of 
primaries. If we associate each primary operator with the certain type of motion then it 
is indeed the simplest (but by no means the only) type of turbulence. Let us look at the 
minimal models of the type (2,2N +1), and assume, that vacuum expectation values are 
absent. In this models we have the set of N primary operators [ip s ], with s — 1, ■ ■ ■ N. 
The fusion rules are: 

foM hU = foU + foks-2] + • • • 

s 3 = min (si + s 2 - 1, 2N + 1 - ( Sl + s 2 - 1)) (68) 
The dimensions are given by: 

(2Af -.)(.-!) 

2(2iV+l) [ ' 

If we identify the stream function ip with the primary field ips with some S , then the flux 
condition gives: 

A s + A 2s _x = -3 (70) 

(if 2 (2s — 1) < 2N + 1). This diophantine equation has unique solution, s=4; N=10 . 
Thus, we conclude that the minimal turbulence is described by the (2 ,21) minimal model. 
Anomalous dimensions in this case are 

A^ = A 4 = - | (71) 
A^ = A 7 = - H (72) 
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As a result we obtain in this case the energy spectrum 

25 

E (k) oo 



(73) 



It seems to be consistent with observations. [7] 

An interesting question if this model is parity conservation. There is no a priori reason 
why we should exclude parity breaking solutions from the consideration. Let us look what 
is the situation in our model. The stream function ^ is a pseudoscalar. That means that if 
parity is conserved, the 3 - point function must contain vector products. But in CFT this 
is impossible: all 3-point functions are simple powers of distances.^ Hence, in conformal 
turbulence the necessary condition for parity conservation is: 

< ^ ^ $ > =0 (74) 

This is indeed the case, since: 

WW-M + W + W + W (75) 

and there is no [i/j^] in the RHS. As we look at higher correlations the situation is less 
clear. Indeed, the 5 - point functions of ip are still zero (by the similar arguments) while 
the 7 -point functions are non-zero. If parity is conserved, they must contain vector 
products. It is not obvious, whether this can be done, or, in other words, whether we can 
combine conformal blocks in the 7 - point function in a antisymmetric (under reflections) 
way. This we leave for the future analyses, noticing only, that if we consider the standard 
symmetric combinations of the conformal blocks, we get a very peculiar picture in this 
model, namely, the parity violation, if present at all, is well hidden - one must go to the 
seven - point functions in order to notice it! 

Let us point out that parity is automaticaly conserved in the non-minimal models. In 
this case there exists infinite number of the degenerate operators, forming closed algebra 
preserving parity [4]. 

7 Conclusion. 

There are infinitely many other solutions of the flux conditions, the fixed points with 
larger number of structures. They have been considered in [8]. It is hard to say, whether 
all matching conditions have been exploited and whether we indeed have infinite number 
of discrete fixed points. If so, the turbulent fluid has infinite discrete regimes and must 
be able to jump from one regime to the other, generating new quasistable structures, we 
are not in a position yet to prove this conjecture, but it is certainly within the scope of 
our methods. 

3 Some time ago A. Migdal conjectured that in 3d turbulence T- invariance forbids the 3-point functions. 
Whether it is true is an interesting open question. In our case the only thing which can be proved is the 
combined PT-invariance . 
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As far as observations are concerned, the best way to check conformal invariance is to 
study momentum representation of the correlation functions when all the momenta and 
all their partial sums belong to the inertial range. In this way we avoid the infrared region 
which is scarcely understood. 

The situation is analogous to that in quantum chromodynamics - there we know 
dynamics at high momenta and quite ignorant in the infrared. Nevertheless, QCD has 
been tested. 

The theory of conformal turbulence gives explicit and unambiguous predictions con- 
cerning these type of correlators. They are expressed in terms of hypergeometric functions, 
and, more importantly have specific "fusion rule" structure when one momentum is get- 
ting much larger than the others. The most direct way to test the theory is to check this 
fusion property. For example, one takes the 4 - point function and reduce it to the 3 - 
point function by taking one large momentum. This 3 - point function in the momentum 
space is a simple hypergeometric function, which becomes just a product of two powers 
if one does the fusion procedure again. This test is especially important, since we are not 
certain which solution of conformal turbulence is realized in a given experiment. Since 
the structure of fusion rules define the theory the test is capable to answer this question. 
On the other hand, may be the most interesting part of the theory is the infrared one, still 
hidden from us. With further work it might become possible to get equations, governing 
the behavior of the vacuum expectation values of < O n > in space and time. This is the 
problem analogous to the one of equation of state in critical phenomena. In any case, the 
theory presented here is but a first step in the long future investigations. 
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